Raman scattering through a metal-insulator transition 
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The exact solution for nonresonant A\ s and Bi g Raman scattering is presented for the simplest 
model that has a correlated metal-insulator transition — the Falicov-Kimball model, by employing 
dynamical mean field theory. In the general case, the Ai g response includes nonresonant, resonant, 
and mixed contributions, the Bi s response includes nonresonant and resonant contributions (we 
prove the Shastry-Shraiman relation for the nonresonant Bi g response) while the i?2g response 
is purely resonant. Three main features are seen in the nonresonant Bi g channel: (i) the rapid 
appearance of low-energy spectral weight at the expense of higher-energy weight; (b) the frequency 
range for this low-energy spectral weight is much larger than the onset temperature, where the 
response first appears; and (iii) the occurrence of an isosbestic point, which is a characteristic 
frequency where the Raman response is independent of temperature for low temperatures. Vertex 
corrections renormalize away all of these anomalous features in the nonresonant Ai s channel. The 
calculated results compare favorably to the Raman response of a number of correlated systems on 
the insulating side of the quantum-critical point (ranging from Kondo insulators, to mixed-valence 
materials, to underdoped high-temperature superconductors). We also show why the nonresonant 
B\ g Raman response is "universal" on the insulating side of the metal-insulator transition. 
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I. INTRODUCTION 

Raman scattering has been an important experimental 
tool for studying lattice dynamics for over four decades 
(since the advent of the laser). More recently, it has been 
applied to study the scattering of electrons in metals, 
insulators, semiconductors, and superconductors. Via 
light's coupling to the electron's charge, inelastic light 
scattering reveals symmetry selective properties of the 
electron dynamics over a wide range of energy scales and 
temperatures. It is similar to the optical conductivity, 
which involves elastic scattering of light by electron-hole 
pairs; but Raman scattering provides additional informa- 
tion, since it is able to isolate different symmetry channels 
by selectively polarizing the incident light and measuring 
the reflected light with a polarized detector. Three prin- 
ciple symmetries are usually examined: (i) A\ g which has 
the full symmetry of the lattice (i.e. is s-like); (ii) 2? lg 
which is a d-like symmetry (that probes the Brillouin- 
zone axes); and (iii) i?2g which is another d-like symmetry 
(that probes the Brillouin-zone diagonals). In 1930-91, 
Shastry and Shraiman proposed a simple relationEI that 
connects the nonresonant Raman response to the optical 
conductivity. We prove the Shastry-Shraiman relation 
here for the B\ g channel in the large-dimensional limit. 

Strongly correlated systems as disparate as mixecL 
valence compoundsp (such as SmBg), Kondo- insulatorsH 
(such as FeSi), and ther-nnderdoped cuprate high temper- 
ature superconductorscl 13, show temperature-dependent 
Big Raman spectra that are both remarkably similar 
and quite anomalous. This "universality" suggests that 
there is a common mechanism governing the electronic 
transport in correlated insulators. As these materials 
are cooled, they all show a pile up of spectral weight for 



moderate photon energy losses with a simultaneous re- 
duction of the low- frequency spectral weight. This spec- 
tral weight transfer is slow at high temperatures and 
then rapidly increases as the temperature is lowered to- 
wards a putative quantum critical point (corresponding 
to a metal-insulator transition). In addition, the Raman 
spectral range is divided into two regions: one where the 
response decreases as T is lowered and one where the 
response increases. These regions are separated by a so- 
called isosbestic point, which is defined to be the charac- 
teristic frequency where the Raman response is indepen- 
dent of temperature. Finally, it is often observed that the 
range of frequency where the Raman response is reduced 
as T is lowered is an order of magnitude (or more) larger 
than the temperature at which the low-frequency spec- 
tral weight disappears. These anomalous features are not 
typically seen in either the Ai g or the £?2g channels. 

Theory has lagged behind experiment for electronic 
Raman scattering in strongly correlated materials. While 
theories that describe Raman scattering in weakly .corre- 
lated (Fermi-liquid) metalsQ or in band insulatorsB have 
been known for some time, it is only recently that a 
theory that describes materials near the metal-insulator 
transition has been developed!! This theoretical treat- 
ment involves applying the dynamical mean-field theory 
to the simplest many-body system that has a quantum- 
critical point — the spinless Falicov-Kimball modeO. We 
choose this model as our canonical model for Raman scat- 
tering because it can be solved exactly and the results are 
universal for the nonresonant B\ g channel on the insulat- 
ing side of the transition (we do not choose it because we 
believe the Falicov-Kimball model is the physically rel- 
evant model for all correlated experimental systems — it 
isn't). Other work has been performed in the so-called 



1 



iterated-perturbation-theory approximation to the Hub- 
bard model, where the nonresonant B\ g Raman response 
is determinedtil. 

The Hamiltonian contains two types of electrons: itin- 
erant band electrons and localized (d or f) electrons. 
The band electrons can hop between nearest neighbors 
[with hopping integral t*/(2Vd) on a d-dimensional cu- 
bic latticeEj], and they interact via a screened Coulomb 
interaction with the localized electrons (that is described 
by an interaction strength U between electrons that are 
located at the same lattice site). All energies are mea- 
sured in units of t* . The Hamiltonian is 

H = --^n E d \ d 3 + E f E w > - »YSA*i + w i) 
+uj24diWi, (i) 

i 

where d\ (di) is the spinless conduction electron creation 
(annihilation) operator at lattice site i and uii = or 1 
is a classical variable corresponding to the localized /- 
electron number at site i. We will adjust both Ef and 
\x so that the average filling of the d-electrons is 1/2 and 
the average filling of the /-electrons is 1/2 (/i = U/2 and 

In this contribution we will show how to derive the non- 
resonant, mixed, and resonant Raman response of a sys- 
tem that crosses through a metal-insulator transition by 
solving the problem exactly in the large dimensional limit 
(employing dynamical mean- field theory). In the case of 
nearest-neighbor hopping on a hypercubic lattice in infi- 
nite dimensions, we show that the A± g response includes 
contributions from all processes, the B\ g response is reso- 
nant or nonresonant, and the £>2 g response is purely reso- 
nant. We provide computational results only for the non- 
resonant Ai g and Bi g responses. Resonant (and mixed) 
Raman scattering results will be presented elsewhere. We 
also prove the Shastry-Shraiman relation, motivate the 
origin of the anomalous features in the nonresonant Bi g 
response, and show why they are not seen in the nonres- 
onant Ai g channel. 

We present the formalism in Section II which includes 
a detailed analysis of the Raman response in the different 
symmetry channels. In Section III we present the results 
of our calculations which display anomalous behavior for 
Bi g and ordinary behavior for Ai g . Our conclusions are 
presented in Section IV. 

II. FORMALISM 

A. Nonresonant Raman response 

The Feynman diagrams for the nonresonant Raman 
response are shown in Fig. [I]. The straight lines denote 
the momentum-dependent Green's function and the wavy 
lines denote the photon propagator. The shading is used 
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FIG. 1. Dyson equation for the nonresonant Raman re- 
sponse function. Solid lines denote electron propagators and 
wavy lines denote photon propagators. The shading denotes 
the fully renormalized susceptibility and the symbol V is the 
irreducible frequency-dependent charge vertex. The vertex 
function 7 is the Raman scattering amplitude, which deter- 
mines the symmetry of the Raman scattering channel. 

to represent the renormalized susceptibility and the sym- 
bol r denotes the local irreducible charge vertex. The 
functions 7(k) are called the Raman scattering ampli- 
tudes, and they are chosen to have a well-defined spatial 
symmetry and no frequency dependence. The Raman re- 
sponse is found from this frequency-dependent density- 
density correlation function which is 

XRaman(*^) = V> / d re WlT 
k J ° 

l Trr r ( e -^p k (r)p k (0)) 

"Tr( e -^p k (0)>l 2 l , . 



with the uniform (q = 0) Raman density operator 

Pk = 7 (k)d k d k , d k = l^ e -^ k ^., (3) 

3 

Z = Tr(e^ l3H ), the partition function, and iv\ = HitlT 
the Bosonic Matsubara frequency (the r-dependence of 
the operators is with respect to the full Hamiltonian). 
The Raman response is characterized in terms of the 
different spatial symmetries of the scattering amplitude. 
One can expand this function in a Fourier series and ex- 
amine the contributions of the lowest components of the 
series, and compare them to experiment. Alternatively, 
one can start from a metallic Hamiltonian, and expand 
the Raman response in powers of the electronic vector 
potential. Under the assumption that the photon wave- 
length is much larger than the lattice spacing, this latter 
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approach yields the following results for the Raman scat- 
tering amplitudes in two dimensions: 



7A lg (k) 
7Bi g (k) 
1b 2s (k) 



d 2 e(k) d 2 e(k) 
dk% dk^ 

d 2 e(k) d 2 e(k) 
dk 2 x 

d 2 e{k) 



dkl 



-<k), 

cos k x — cos k y , 



sin fc x sin k„ 



(4) 



cos fey the electronic bandstructure. 



with e(k) cx cos k x 
Note that the B2 g response vanishes for pure nearest- 
neighbor hopping, which is what we consider here (the 
band structure is just a sum of cosines; hence all mixed 
derivatives vanish). The nonresonant B 2g Raman re- 
sponse then vanishes, because the lowest-order Raman 
scattering amplitude vanishes. 

The above forms can be generalized to the infinite- 
d limit -for non-resonant scattering by noting that 7(k) 
satisfies^ 



7(k) = E' 



a/3 



8 2 e(k) 

dk a dkp 



(5) 



with (e c ) the incoming (outgoing) photon polarization. 
In systems with hopping beyond nearest neighbors, it is 
not possible to project completely onto the A\ & sector 
with just one choice of polarizations, since either B\ g 
or £?2g sectors will always be mixed in; it is possible to 
project onto i?i g or B 2g with one measurement. In the 
infinite-dimensional limit, with nearest-neighbor hopping 
only, we can choose ei a = e oa = 1 for Ai g , ej Q = 1, 

G oa ( 1) for -Big, and G^q— even 0; ^ia — odd 1; 

e OQ=cvcn = 1, e OQ=oc jd = for B 2g . The resulting Raman 
scattering amplitudes are 

d 

7Alg (k)«c-e(k), 7Slg (k)«-=]T(-ircosfc 4 , (6) 

^ d i= i 

and (k) = 0. Note that we include a constant term 
c in the Ai g amplitude, since it is allowed by symmetry. 

The nonresonant Ramap_Eesponse in the B lg channel 
has no vertex correctionsOEl and is equal to the bare 
bubble. The key element needed to show this is that the 
irreducible charge vertex T is independent of momentum 
(since it is local). Hence, we must evaluate a summation 
over k of the form 
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d 



_ , (— lVcos/ci 



X 



(7) 



which arises when solving the Dyson equation for the 
Big response (and we can assume without loss of gen- 
erality that the imaginary part of X is greater than 
zero). Writing the fraction in Eq. (R) as the integral 
of an exponentiate 



X 



2 v^d 



cos ki 



poo 2 

'. / dzexp[iz(X-\ — —\ cos kj)], 
Jo 1 v^^t 



(8) 



allows us to decouple the summation over momentum to 
the sum over d identical terms, each multiplied by (— 1)\ 
This then vanishes for all even d and for odd-d in the limit 
d — * oo (due to the I /yd term). So the nonresonant B\ s 
response reduces to the evaluation of the bare bubble. 

In the j4i g channel, the Raman scattering amplitude 
has the same symmetry as the irreducible charge vertex, 
so the corresponding summation over momentum does 
not vanish, and the nonresonant A\ g Raman response is 
renormalized by the irreducible charge vertex. 

The Falicov-Kimball model can be solved exactly in 
the infinite-dimensional limit by using dynamical mean- 
field theory (see Ref. [l6] for details). We summarize the 
main points to establish our notation. The local Green's 
function at the Fermionic Matsubara frequency iu> n — 
inT(2n + 1) is defined by 



G n = G{iu! n ) 



-0H a 



-TrT T / dre 1 ^ 
'o 



d(T)d^{0)S(X)) 



Z 



(9) 



with 



Z = Z (ti) + e-^ E f-^Z ( t i-U), (10) 

the atomic partition function expressed in terms of 

Z (fi) = Tr d {e-P Ha S(X)), H = -fidU. (11) 

In the above equations, the atomic Hamiltonian H at is 
the Hamiltonian of Eq. (|l|) restricted to one site, with 
t* = 0, and all time dependence is with respect to H at . 
The evolution operator 5(A) satisfies 

rP rP 

S{\) = cxp[- dr dr' ' d) (t)A(t - r')d(r')], (12) 
Jo Jo 

with A(r — t') a time-dependent atomic field adjusted 
to make the atomic Green's function equal to the local 
lattice Green's function. We define an effective medium 
by 



G 1 (iu n ) = G n 1 + £„ = iu n +n- A r 



(13) 



with S ra the local self-energy and X n the Fourier trans- 
formation of A(t). The trace in Eq. (||) can be evaluated 
directly to yield 

G n = waGo{iuJn) + w^Gq 1 (iu„) - U}~\ (14) 

with Wq = 1 — Wi, 

w 1 =exp[-p(E f -ij)]Zo{iJL-U)/Z, (15) 
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and 



i \ 

/2 TT lLU n + M — "V. 



(16) 



n— — oo 



The self-consistency relation needed to determine A„ and 
G n is to equate the local lattice Green's function to the 
atomic Green's function via 



G n — 



de- 



Pie) 



lOJ n 



(17) 



with p(e) = exp(— e 2 )/y / 7r the noninteracting density of 
states for the infinite-dimensional hypercubic lattice. 

The iterative algorithm to solve for G n starts with 
E n = 0. Then Eq. is used to find G„, Eq. (|1) 
is employed to extract the effective medium, Eq. (14) 
is used to find a new local Green's function, and then 
Eq. (|r]) is used to find the new self-energy. This se- 
quence of steps is then repeated until it converges, which 
usually requires only about a dozen or so iterations. This 
algorithm can also be used on the real axis (with suitably 
modified equations) to directly solve for the Green's func- 
tion and self-energy on the real axis. 

Once the Green's functions are known, then the non- 
resonant Raman response can be calculated directly. The 
Bi g channel is simple, since we need only evaluate the 
bare bubble. Noting that the average of the square of 
the Bi g Raman scattering amplitude is 1/2 then yields 



XB ls (m) 



T 



E 



G(iui n ) - G(iu n+ i) 
in + E(iu>„) - T>(iu> n+ i) ' 



(18) 



for the Raman response on the imaginary axis. This for- 
mula can be easily analytically continued to the real axis 
by following the same procedure outlined in—the calcu- 
lation of the dynamical charge susceptibilitycll: rewrite 
the sum over Matsubara frequencies by a contour inte- 
gral of advanced or retarded Green's functions and self- 
energies multiplied by the Fermi factor, and then deform 
the contours to the real axis picking up any poles in the 
complex plane. Under the assumption that there are no 
extra poles when the contours are deformed, one ends up 
with the following expression for the Bi g response: 



XBigO) 



4tt 



du> 



G(u) - G(lu + v) 



f(u + u) 



v + E(w) — E(k> + v) 

G* (lu) ~ G* (uj + v) 



v + E* (lu) — E* (lu + v) 



-[/m -/(„+„)] \, d9) 

v + E*(w) — Zj(uj + v) I 

with f(u>) = l/[l + exp(/3o;)] the Fermi function. We ver- 
ify that this expression is indeed accurate, by using the 
spectral formula to calculate the Raman response on the 



imaginary axis and comparing it to the result directly cal- 
culated from the expression in Eq. (p!s|) . We find that the 
results rarely differ by more than one part in a thousand 
confirming the accuracy of the analytic continuation (we 
believe, but cannot prove, that no additional poles ex- 
ist in the complex plane, which would mean the analytic 
continuation is exact). 

We can use the expression in Eq. (|l9| ) to prove the 
Shastry-Shraiman relation. Using the identity that re- 
lates the local Green's function on the lattice to the self 
energy in Eq. ( |l7j ) (but evaluated on the real axis rather 
than the imaginary axis) yields 



G(lo) - G(lo + v) 
v + E(w) — E(a; + v) 



de-. 



+ fj, — E(cij) — e] [lu + v + fi — E(oj + v) — 



(20) 



Employing this identity, and related ones for the ad- 
vanced Green's functions, in Eq. (|1]) produces 

XB ls (v) = -zz I duj dep(e) 

^ J —oo J —oo 

f m 



Im- 



LO + v + [i — E(a; + v) — e lu + \i — E(w) — e 



f(LU + V) 



-Im 



lu + fi — E*(w) — e lu + v + /i — Y,(lu -\- v) — e 
Using the definition of the spectral function 
A(e,Lu) = - -Im- ' 



(21) 



7T LU + fi — E(cij) — 6 ' 

and taking the imaginary part of Eq. ( pl| ) , then yields 



(22) 



/oo poo 
dcu / dep(e) 
-oo J — oo 



[ f(u)-f(u + v)]A(e,u)A(e,u + v), 
cx va(v), (23) 

which is the Shastry-Shraiman relation — the imaginary 
part of the nonresonant i?i g Raman response is propor- 
tional to the frequency v times the optical conductivity 
o(v). This conclusion comes from comparing the integra- 
tion in Eq. ( p3| ) to the well-known result for the-pptical 
conductivity on the infinite-d hypercubic latticeEa. Note 
that this final formula for the Raman response depends 
only on the shape of the spectral function. Since the 
derivation was modcl-indppcndcnt (like the derivation of 
the optical conductivityca) , this form for the nonresonant 
Big Raman response holds for all models. (Similarly, in 
models with next-nearest-neighbor hopping, one can also 
show that the nonresonant B2 g response also satisfies the 
Shastry-Shraiman relation.) We note that this is a conse- 
quence of the momentum-independent self energy — any 
dependence of the self energy on momentum will gener- 
ally violate the Shastry-Shraiman relation, although it 
may still be an accurate approximation. 
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The Ai g response is more complicated, because it re- 
quires a proper treatment of the vertex contributions. 
Fortunately, the chacge vertex for the Falicov-Kimball 
model is well-knowr£3 and assumes a simple form (for 

m + o) 



T(iu> m ,iu) n ;ii'i) 



1 Y,(iu n ) - T,(iuj n+ i) 
l TG(iw n ) - G(iuj n+ i)' 



(24) 



Hence, the Raman response in the A lg channel can be 
found by solving the relevant Dyson's equation, using 
the above form of the charge vertex. The set of Feynman 
diagrams is shown in Figs. [I] and ||. Note that we have 
to solve these two coupled equations for the A\ g Raman 
response. 




- T 
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FIG. 2. Supplemental Feynman diagrams for the nonreso- 
nant A lg Raman response. These diagrams are identical to 
those in Figure 1 except they have one fewer power of 7. 

The difference between the two diagrams shown in 
Figs, [l] and || is the number of factors of the Raman 
scattering amplitude 7; the nonresonant Raman response 
requires the dressed response that has only one power of 
7 — this dressed response satisfies a simple Dyson equa- 
tion, shown in Fig. ||. We denote the dressed response 
function in Fig. || by x' {iu m , iu n ]iv{). Then a straightfor- 
ward evaluation of the Feynman diagrams and a solution 
of the Dyson equation produces 



'-''in/a 



1 + xo{i^ m ;m)TT(iLu mi iuj m ;ivi)' 



(25) 



with the irreducible charge vertex function found in 
Eq. (pi), 



xo(iu n ;m) = - 

and 



G(iuj n ) - G(iuj n+ i) 



ivi + T,{ioj n ) - Yj(iijj n+ i) 1 

c(G n — Gn+l) ~ Z n G n + Zn+lGn+l 



(26) 



(27) 



Here we used the notation Z n = iu> n + [i — Y,(ioj n ). Now, 
knowledge of x' allows us to solve for the Raman response 
in Fig. |l|. The end result is 

n 

Xo(iLLi n ; ivi) - [xo(w„; iv{) + G n G n+ i]TT(im n , iuj n ; ivi) 
l + Xo{iu n ]ivi)TT{iLu ni iu} n ]ivi) 

where the charge vertex is found in Eq. (p4|), the bare 
susceptibility xo is found in Eq. (|2^) and the other bare 
susceptibility Xo (which is where all of the c dependence 
lies) is 



(28) 



Xa{iu n ;in) = [-c 2 (G n - G n+t ) 

+2c(Z n G n — Z n+ iG n +i)^ 
+ Z n — Z n +i — Z n G n 
/[in + s n - £„+/]. 



(29) 



Surprisingly, all of the c dependence actually drops out 
of the nonresonant Raman response (as it does in 
conventional metalslZl). This can be seen by substituting 
the representations for xo and T into the denominator in 
the first term of Eq. d2q) . One finds that the summation 
for the first term can then be written in the general form 
^nO^n — Y n+ i)/ivi, which vanishes since the individual 
summations converge and one can change the summation 
index n — > n + I in the first term [we similarly note that 
^ n (I] n — £„+;) = 0]. Hence we arrive at the final form 
for the nonresonant A\ g Raman response: 

/ • \ T (S„ — Yi n+ i)G n G n+ i 
XA ls {n>i) = — > T . r-r 

_ T ^--^ (T, n — Y, n+ i){ivi + £„ — S„ + () 



ay 1 



It is a straightforward exercise to perform a similar ana- 
lytic continuation of this expression. The result is 



2-kv 



[S(w) - £(w + v)][v + - S(w + v)\ 



G-^uj + ^-G-^lo) 
[£* (w) - S* (oj + v)] [v + S* (uj) - E* (u + v)] 



-[f(w)-f{w + v)] 

[S*H - Z{LU + !/)][!/ + E*(cj) - S^> + v)] ] 

G-V + i^-G- 1 *^) I 



(31) 



B. Resonant Raman response 

Resonant Raman scattering arises from a fourth-order 
process with respect to the photon vector potential. One 
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(a) (b) 
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FIG. 3. Bare resonant Raman scattering diagrams. The 
(a) uncrossed, (b) first partially crossed, (c) second partially 
crossed, and (d) fully crossed diagrams are all shown. The 
straight lines are momentum dependent Green's functions, the 
wiggly lines are photon propagators, and the vertex factors 
are the dot product of the photon polarization with the Fermi 
velocity as described in the text (vi = ej • v, v = e • v). Note 
that the two partially crossed diagrams (b) and (c) are equal 
to each other. 

way to view this scattering is to "separate" the two- 
photon-two-electron vertex of nonresonant Raman scat- 
tering in Fig. |l| into two single-photon- two-electron ver- 
tices. Keeping in mind the direct and exchange possibil- 
ities, results in four distinct bare resonant Raman scat- 
tering diagrams, which are depicted in Fig. [| (a) the 
uncrossed (direct) diagram; (b) the first partially-crossed 
diagram; (c) the second partially crossed diagram; and 
(d) the fully crossed (exchange) diagram. In all relevant 
cases, the partially crossed diagrams (b) and (c) are equal 
to each other. 

In general, the Raman vertices for resonant scatter- 
ing involve expectation values of the momentum opera- 
tor between the conduction band and the excited states. 
However, in the single-band model, the Raman vertices 
in Fig. H are equal to the dot product of the photon po- 
larization vector e with the Fermi velocity v 

v a = = 2sinfc Q . (32) 

Ok a 

Again, we distinguish between the incoming photon po- 
larization e, and the outgoing photon polarization e Q . 
The different symmetry channels are projected by ap- 
propriate choice of the polarization vectors. We choose 
eia = e OQ = 1 for A\ g , e ia = 1, e oa = (-1)" for B lg , and 



e^a— oven — t), €iia— odd — 1; e oa — cvcn 1, G oa — dd — O tor 

£?2g. It is important to note that in all cases, a single 
factor of e • v is an odd function of k, so the only way to 
get nonzero summations over momentum is to have an 
even number of e • v factors in any given integration. 

In general, the bare resonant Raman scattering dia- 
grams must be renormalized by attaching appropriate ir- 
reducible charge vertex functions (and higher-order gen- 
eralizations, if possible) to all relevant Green's function 
legs. Since all e • v factors are odd in k and the irre- 
ducible charge vertex has A\ g symmetry, there can be no 
renormalization of any single vertex. Similarly, one can 
argue that there are no three-particle or four-particle ver- 
tex renormalizations possible either. The only possibility 
is a two-particle vertex that connects opposite Green's 
function lines. In the Ai g resonant Raman scattering 
channel, all possible vertical and horizontal renormaliza- 
tions of each bare diagram is possible. For the B\ g and 
£?2g channels, there are simplifications. One can quickly 
verify that in both of those cases the product e$ • ve D • v is 
orthogonal to the A\ g symmetry, so the partially crossed 
diagrams (b) and (c) cannot have any renormalization 
either. Furthermore, this also implies that attaching a 
vertical charge vertex to either the (a) or the (d) dia- 
gram vanishes for the same reason. But a horizontal 
attachment of the vertex is possible in both the (a) and 
(d) diagrams. Hence, the renormalized resonant Raman 
scattering in both the Bi g and B^ g channels is described 
by the coupled set of equations in Fig. ^|. A similar, but 
more complicated, set of diagrams is needed for the A\ g 
channel, where renormalizations are present on both the 
vertical and horizontal pairs of legs for all diagrams. We 
don't show those Feynman diagrams here. 

Evaluation of these diagrams, and their analytic con- 
tinuation to the real axis is tedious. We leave that task 
for a future publication. But we do note that both of the 
renormalized diagrams (a) and (d) of Fig. | will have the 
zero-frequency (ivi = 0) piece of the irreducible charge 
vertex renormalizing them. This is the nondiagonal piece 
of the charge vertex, and it can get large when the sys- 
tem is tuned to lie near a phase-separation transition. 
Hence, we expect the resonant Raman scattering to be 
enhanced whenever one is close to a phase separation. 
Note that this enhancement will occur for all photon fre- 
quencies, since the zero-frequency vertex couples all such 
diagrams; this implies that one will not see this effect by 
tuning through a resonant frequency. 

C. Mixed Raman response 

The mixed Raman response comes from the cross terms 
between the linear and quadratic terms in the vector po- 
tential. As such, these diagrams have two single-photon- 
two-electron vertices and one two-photon-two-electron 
vertex. Since the £?2 g Raman scattering amplitude van- 
ishes for nearest-neighbor hopping on a hypercubic lat- 
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FIG. 4. Renormalized resonant Raman scattering diagrams 
for Figures 3(a) and 3(d) in the Bi g and B2 g channels. We 
suppress the photon lines that are explicitly shown in Fig. ^| 
for simplicity. In (a) one has the diagrams for the Raman re- 
sponse, while in (b) the supplemental Dyson equation needed 
to solve for the dressed Raman response is given. The only dif- 
ference for the Big and B2 g channels is the choice of polariza- 
tion vectors e* and e , and we have suppressed the incoming 
and outgoing indices on the vertices v. The symbol F denotes 
the local irreducible charge vertex. Note that only the direct 
and exchange diagrams [Figs. |§](a) and (d)] are renormalized. 
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FIG. 5. Bare mixed Raman response. The single-photon 
vertices are multiplied by the respective factor of the po- 
larization vector dotted into the Fermi velocity, while the 
two-photon vertex is multiplied by the corresponding Raman 
scattering amplitude. Note that no renormalization is possi- 
ble for the Big channel, but the diagram is renormalized in 
the Aig channel. 



tice, there is no mixed Raman response for that chan- 
nel. The bare mixed Raman response is shown in Fig. |^. 
There are two possible diagrams corresponding to the 
direct or exchange processes. In the B\ g channel, the 
dressed mixed response is equal to the bare mixed re- 
sponse, since the irreducible charge vertex cannot be in- 
serted anywhere. For the A\ g channel, the bare mixed 
response is renormalized by the irreducible charge vertex 
in a similar way to how it renormalized the nonresonant 
Raman response. We don't present any numerical results 
for the mixed Raman response here, but will do so in an- 
other publication. It turns out, that one can show that 
in the B\ g channel with nearest-neighbor hopping only, 
the mixed response is a 1/d correction and can be ne- 
glected; it cannot be neglected for the Ai g sector. The 
1 /d correction for B\ g arises from the fact that the sum- 
mation over momentum has a "form factor" proportional 
to cos(fci)(— iy sin 2 (fcj)(— l) 3 . This term cancels when 
summed over i and j except for i = j. This latter con- 
straint forces the mixed diagram to be a 1/d correction. 
In the Ai g case, there are no factors of (—1) so the terms 
with all i and j contribute, and it is an O(l) term. 



III. RESULTS 

The Falicov-Kimball model has a ground state that is 
not a Fermi liquid because the lifetime of a quasiparticle 
is finite at the Fermi energy. In addition, the imaginary 
part of the self energy has the wrong sign of curvature to 
be a Fermi liquid. We study the model at half filling. As 
U increases, the system first enters a pseudogap phase, 
where spectral weight is depleted near the chemical po- 
tential, and then undergoes a metal-insulator transition 
(the pseudogap phase is possible because the ground state 
is not a Fermi liquid). The interacting density of states 
(DOS) is, however, temperature-independent for fixed U 
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FIG. 6. Interacting density of states for the Fali- 
cov-Kimball model. Results are shown for U = 0.5, 1, 1.5, 
2, and 4 (the numbers in the figure label the value of U). 
Note how the system first develops a pseudogap (1.0) before 
the metal-insulator transition at U = 1.5. The density of 
states is independent of temperature. 

and fixed electron fillings^. It is plotted in Fig. ^| for a 
range of values of U: U < 0.65 corresponds to a weakly- 
correlated metal, while a pseudogap phase appears for 
0.65 < U < 1.5 moving through a quantum critical point 
at U = 1.5 to the insulator phase U > 1.5 (we neglect all 
possible charge-density- wave phases here) . 

In Figure CTa) we plot the nonresonant B\ & Raman re- 
sponse at a fixed temperature T = 0.5 for different values 
of U. For small values of U, a small scattering intensity is 
observed due to the weak interaction among "quasipar- 
ticles" providing a small region of phase space allowable 
for pair scattering. The peak of the response reflects the 
dominant energy scale for scattering, as is well known in 
metalso and the high-energy tail is the cutoff determined 
by the finite energy band. This shape is also understand- 
able from the Shastry-Shraiman relation — since the op- 
tical conductivity is a Lorentzian, the Raman response 
is just proportional to av/[v 2 + a 2 ], which assumes the 
above form. As U increases, the low-frequency response 
is depleted as spectral weight gets shifted into a large 
charge transfer peak at a frequency ~ U. The charge 
transfer peak begins to appear for values of U for which 
the DOS is still finite at the Fermi level (U = 1) and be- 
comes large in this pseudogap phase before growing even 
larger in the insulating phase. Notice how low-frequency 
spectral weight remains even as one is well on the insu- 
lating side of the quantum critical point (U — 4) and at 
a temperature T much lower than the gap. It is these 
spectral features that are characteristically seen in the 
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FIG. 7. Nonresonant (a) Bi g and (b) Ai s Raman response 
for different values of U at T = 0.5. The Raman response is 
measured in arbitrary units. Notice how the vertex correc- 
tions suppress the low-frequency spectral weight in the Ai g 
channel. 

experiments and which can only be seen in a theory that 
approaches the quantum critical point. 

In Fig. 0(b), we show the nonresonant Ai g Raman re- 
sponse. As shown above, the A\ g response is independent 
of the value of c. The general behavior is similar to that 
of the Big channel except (i) at weak coupling the Ra- 
man scattering is more symmetric and pushed to higher 
energy; (ii) the vertex corrections suppress all nontrivial 
low-frequency Raman response; and (iii) the widths of 
the charge-transfer peaks are enhanced. 

Since the Raman response displays anomalous features 
on the insulating side of the metal-insulator transition, 
we first present results for U — 2, just on the insulating 
side of the quantum critical point. In Figure ^(a) , we plot 
the temperature dependence. The total spectral weight 
increases dramatically with decreasing temperature as 
charge transfer processes become more sharply defined. 
At the same time, the low-frequency response depletes 
with lowering temperatures, vanishing at a temperature 
which is on the order of the T — insulating gap (we 
are unable to analytically estimate the crossover temper- 
ature). This behavior is-precisely what is seen in exper- 
iments ontl FeSi and onQ underdoped La2_ a; Sr :c Cu04 at 
low temperatures where both the isosbestic point and the 
low temperature spectral weight depletion can be seen. 
Similar results are also seenH in SmBg, but a low-energy 
peak also develops in that material at low temperatures. 
In the Ai g channel [Fig. ||(b)] , we see the same sharpening 
of the charge-transfer peak at low T, but the low-energy 
response is much smaller and changes much more slowly 
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FIG. 8. (a) Nonresonant Bi g Raman response for a range 
of temperatures (T = 0.05, 0.2, 0.3, 0.5, 0.9) for U = 2 (which 
lies just on the insulating side of the metal-insulator transi- 
tion) and (b) nonresonant A\ % Raman response for the same 
temperatures at U = 2. The lines are labeled by their tem- 
perature (except for T = 0.2 which is unlabeled). Note how 
the Big response has low-frequency spectral weight that de- 
velops rapidly at an onset temperature of T « 0.2 (the low 
frequency response at T=0.5 and T=0.9 overlap) and note the 
isosbestic point at v w 1. The ratio of the range in frequency 
over which the low-frequency weight increases and the onset 
temperature is about 5. There are no anomalous features in 
the Ai g spectrum. 



with T (but, in fact, increases as T is lowered). 

If one were to interpret the temperature at which the 
Big Raman spectral weight starts to deplete as the "tran- 
sition temperature" T c and the range of frequency over 
which the weight is depleted as the gap A, then one 
would conclude that near the quantum critical point 
2A/k B T c > 1. This is because the "T c " is effectively 
determined by the gap in the single-particle density of 
states (which is small near the quantum critical point), 
while the "A" is determined by the width of the lower 
Hubbard band (which remains finite at the quantum crit- 
ical point); hence the ratio can become very large near 
the quantum critical point (and should decrease in the 
large- U limit). 

The nonresonant Raman response is plotted in Fig. ^ 
for a number of different temperatures at U = 1. Note 
how the Bi g response has nontrivial low-energy spec- 
tral weight, even though it is not completely separated 
from the charge-transfer processes. Even in this case, one 
can see the low-temperature development of an isosbestic 
point near U/2 for temperatures below about T = 0.3. 
As the low-energy spectral weight is depleted, the peak 
becomes more symmetric in shape. In the A\ & chan- 
nel, the response sharpens, and the peak moves to lower 
energy as the temperature is lowered. In fact, the low- 
energy spectral weight actually increases as T is lowered. 
There is no indication of a separation of the response 
into low and high-energy features that have a different 
temperature dependence (as seen for the B\ g response). 

The Big spectral-weight transfer from low frequencies 
to the charge transfer peak as a function of tempera- 
ture can be quantified by separating the Raman response 
into two regions determined by the isosbestic point and 
plotting the ratio of the low-frequency spectral weight at 
temperature T to the low-frequency spectral weight at 
T = 0.95 versus reduced temperature T/0.95 [Fig. [To| |. 
Choosing U/2 as the location of the isosbestic point 
(which divides the low-frequency and high-frequency re- 
gions), we find that the reduction of spectral weight from 
high to low temperatures is over 50 percent even in the 
weak pseudogap phase, and decreases by well over three 
orders of magnitude as U increases into the insulating 
phase (U = 4). 

The origin of the isosbestic point in the nonresonant 
Big response is mysterious, but can be motivated by the 
Shastry-Shraiman relationEI. Since the optical conductiv- 
ity satisfies a sum rule, the appearance of an isosbestic 
point there is not surprising, as any decrease in low fre- 
quency spectral weight must be compensated by a cor- 
responding increase in high-frequency spectral weight, so 
one might expect there to be a frequency where the re- 
sponse doesn't depend on T (of course this doesn't estab- 
lish that an isosbestic point must exist, it just motivates 
such an existence). The isosbestic point in the Bi g Ra- 
man response then follows from the Shastry-Shraiman 
relation, since multiplying the optical conductivity by a 
frequency will not modify the appearance of an isosbestic 
point. 
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FIG. 9. Nonresonant (a) Bi g and (b) yli g Raman response 
at (7 = 1 for T =0.9, 0.5, 0.3, 0.1, 0.05. The lines are labeled 
by the temperature. Note how the B\ s response develops an 
isosbestic point at low temperatures, but the low-frequency 
depletion is more modest here than in the insulating phase. 
The A\ s response has no anomalous features, but the peak of 
the response moves to lower energies as the Raman response 
sharpens at low temperature. In the moderate-to-low fre- 
quency range, the Raman response increases as T is lowered. 
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FIG. 10. Ratio of low-frequency spectral weight at tem- 
perature T to the low-frequency spectral weight at T = 0.95 
plotted versus reduced temperature T/0.95 for different val- 
ues of U. The values of U are 0.5,1,1.5,2, and 4, and the 
curves correspond to in creasing values of U starting at the 
top with 0.5 and running to the bottom at 4. Note how the 
sharpening of the Raman response as T is lowered results 
in significant reductions to the low-frequency spectral weight 
even in the metallic case. In the strong insulator phase, the 
spectral weight decreases by over three orders of magnitude 
((7 = 4). 



We attribute the presence of anomalous low-frequency 
and low-temperature response in a system which is a 
strongly correlated insulator to the appearance of ther- 
mally activated transport channels (indeed, the only tem- 
perature dependence to the Raman response comes from 
Fermi factors in an integral). In the insulating phase at 
zero temperature, the only available intermediate states 
created by the light must involve double site occupancy of 
a conduction and a localized electron, with an energy cost 
of 7. This gives the large charge transfer peak. As the 
temperature is increased, double occupancy can occur 
and as a result light can scatter electrons to hop between 
adjacent unoccupied states either directly or via virtual 
double occupancies. The number of electrons which can 
scatter in this fashion increases with increasing tempera- 
ture, leading to an increase in the low-frequency spectral 
weight. The frequency range for this low- frequency Ra- 
man response is determined by the lower Hubbard band- 
width, which is typically much larger than the tempera- 
ture at which these features first appear. 

There is also a more mathematical explanation to the 
low-frequency spectral response. If we examine the in- 
tegral for the Big Raman response in Eq. (|l9|), we note 
three important points (i) the imaginary part of the Ra- 
man response is proportional to the real part of the inte- 
grand [the terms within the braces in Eq. (fl9)]; (ii) the 
integrand vanishes if the Green function (and self-energy) 
at to and uj + v are both real; (iii) all temperature depen- 
dence arises from the Fermi factors, since both G and £ 
are temperature-independent. In the insulating regime, 
the DOS breaks into two pieces, a lower subband centered 
at —7/2 with a width of 0(1) and an upper subband at 
7/2 with a width of O(l). The Green's functions are 
complex only when the frequency argument lies within 
one of the bands. Hence there are two main contribu- 
tions to the Raman response: (i) intraband processes, 
where u> w —7/2 or 7/2 and v w 1; and (ii) interband 
processes, where u> ~ —7/2 and v w U. The interband 
processes, with v sa 7 are what give rise to the charge- 
transfer peaks seen in the Raman response; these pro- 
cesses survive even at T = 0. The intraband processes, 
with v 1=3 1, give rise to the low-frequency spectral fea- 
tures. At low temperatures, these features are propor- 
tional to j(to) — f(uj + v) which can be approximated 
by exp(— 7/2T)[l — exp{y/T)] in the insulating phase. 
Hence, we expect the low-frequency spectral weight to 
be proportional to exp(— U/2T) in the large-7 limit. In 
the Ai g channel, the charge vertex makes the integrand 
more complicated to analyze (but the response still sep- 
arates into interband and intraband processes), and the 
vertex corrections end up ultimately suppressing the low- 
frequency response. 

Lastly, we plot the inverse slope of the Raman response 
in Fig. |llj as a function of temperature for different values 
of 7. The inverse Raman slope is the reciprocal of the 
slope of the Raman response in the limit as u> — > 0. Since 
the self energy is temperature independent, we might 
expect a constant Raman slope as a function of tem- 
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FIG. 11. Inverse Raman slope as a function of temperature. 
The inverse Raman slope measures the scattering rate of the 
"quasiparticles" of a correlated metal. As the system enters 
the pseudogap phase, the inverse slope starts to increase at 
low temperatures, increasing dramatically as one enters the 
correlated insulator (U = 2). 

perature, as is the case with a disordered Fermi liquid. 
However, this is not the case due to the formation of a 
thermally generated band for scattering. For small values 
of U, the temperature dependence of the Raman inverse 
slope is weak due to the temperature independence of 
the self energy. However, as the single-particle bands 
begin to separate, the relevance of thermally generated 
double occupancies becomes more pronounced and the 
inverse slope rapidly rises at low temperatures. As the 
system becomes more insulating, the low temperature in- 
verse slope increases dramatically due to the depletion of 
the low- frequency spectral weight. As U increases from 
the pseudogap phase into the insulating phase, the tem- 
perature dependence of the Raman inverse slope indi- 
cates the formation of gapped excitations (assuming the 
form T[l + cosh(A/{2T})], with A the gap in the single- 
particle DOS). Such behavior has been seen in the un- 
derdoped cuprate materials. 



IV. CONCLUSIONS 

The electronic Raman response of a wide variety of 
correlated materials (on the insulating side of the metal- 
insulator transition) display similar anomalous features, 
which point to a common explanation. In our results, 
we have shown how to see these anomalies in the B\ s 
channel, by solving the Falicov-Kimball model. We saw 
that the Raman response ultimately depended only on 
the single-particle density of states. In more complicated 
correlated models, the metallic single-particle density of 



states can have an additional Fermi liquid peak at low 
frequencies, which will add new features to the Raman 



response,EJ but on the insulating side of the transition, 
where most of the anomalous behavior is seen, the single- 
particle density of states must be similar to that of the 
Falicov-Kimball model (except for some additional weak 
temperature dependence of the interacting DOS), since 
an insulator has no low-energy spectral weight. Hence, 
these anomalous Raman scattering results are expected 
to be essentially model-independent (since they only de- 
pend on the interacting DOS)! 

In this work, we also illustrated how one can calculate 
both the resonant and the mixed Raman scattering re- 
sponse as well. We showed how the bare diagrams are 
renormalized for the different symmetry sectors, but did 
not perform any numerical calculations here. 

Our theoretical results compare quite favorably to the 
experimental results seen in a variety of different mate- 
rials rangiiig from mixed-valence compoundso, to Kondo 
insulatorsa to thc^u-odcrdoped high-temperature super- 
conducting oxidesuul. The experimental data illustrate 
the three characteristic features seen in our theory: (i) 
there is a rapid rise in the low-frequency spectral weight 
at low temperatures (at the expense of the high-frequency 
spectral weight), (ii) there is an isosbestic point, and 
(iii) the range of frequency over which the low-frequency 
weight appears is much larger than the onset tempera- 
ture, where it is first seen. Our model always produces 
an isotropic gap, so we are unable to illustrate the sym- 
metry selective behavior seen in the copper-oxides where 
only the B\ g response is anomalous, and the A\ g and 
Big responses are metallic. But our results do indicate 
a "universality" and model independence of the Raman 
response on the insulating side of, but in close proximity 
to, a quantum critical point. We believe this is the reason 
why so many different materials show the same generic 
behavior in their electronic Raman scattering. 

In future work we will examine nonxesonant B\ g Ra- 
man scattering in the Hubbard modelEJ and will exam- 
ine resonant and mixed Raman scattering effects in the 
Falicov-Kimball model. 
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